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00 I Abstract 

Recent research shows that the partition function for a class of 
^ I models involving fermions can be written as a statistical mechanics 

. of clusters with positive definite weights. This new representation 

. of the model allows one to construct novel algorithms. We illus- 

I frate this through models consisting of fermions with and without 

spin. A Hubbard type model with both attractive and repulsive 
. interactions becomes tractable using the new approach. Precision 

results in the two dimensionai attractive model confirm a superfluid 
phase transition in the Kosterlitz-Thouless universality class. 



Oh! 

r-| ! 1 Introduction 



Fermion algorithms are known to be notoriously difìicult. The main reason 
for this is that it is difìicult to write the partition function of models in- 
5^ I volving fermions as a sum over confìgurations with positive definite weight. 

The most common approach is to integrate out the fermions in favor of a 
fermion determinant which leads to an effective bosonic partition function. 
Typically one obtains 

Z = J [d(j)] exp{-S[(f)]) Det(M[(/)]) (1) 

where the determinant is a non-local function of the bosonic fìelds [0] . In 
cases where Det(M[0]) is positive useful algorithms can be found [1-3]. In 
other cases one can make progress only by using uncontrolled approxima- 
tions 1^. Thus it is important to fìnd alternative approaches to fermionic 
path integrals. 
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Recently a novel approach to solve certain lattice fermionic models was 
discovered [5-11]. It is possible to rewrite the partition function of these 
models as a sum over confìgurations of locai bond variables with positive 
definite weights. The bonds b connect lattice sites and thus divide the 
lattice into clusters. The partition function can be written as 

Z ^ J2W[b]Si^[b] (2) 

[b] 

where W[b] > is the magnitude of the Boltzmann weight for a given 
bond configuration, and Sign[ò] > is an entropy factor arising from 
"cluster flips" that result due to degrees of freedom other than the bonds" . 
This representation of the fermionic partition function allows one to con- 
struct cfficient cluster algorithms which had prcviously been found only 
for bosonic problems. In this article we describe this new approach to 
fermionic path integrals and present a result from recent Monte-Carlo 
studies using the new method. 



2 Fermion World-Line Path Integrals 

Consider spin-less fermions, hopping on a rf-dimensional cubie lattice con- 
sisting oì V = L'^ sites and satisfying periodic or anti-periodic spatial 
boundary conditions. Let us focus on models whose Hamilton operators 
can be described by 

H^Y.K,.^i (3) 

x,i 

where i — 1,2, ..d represents directions and hx^i is a nearest neighbor 
operator made up of the usuai fermion creation and annihilation operators 
c+ and Cx associated with the site x. In order to write a path integrai for 
such a problem, the Hamilton operator is decomposed into 2d terms 

H^H^ + H2 + ... + H2d, (4) 

with 

x=(xi,X2,...,Xjj) x=(xi,X2,...,a;j;) 

Xi even odd 

Note that the individuai contributions to a given Hi commute with cach 
other, but two different iJj's do not commute. Using the Suzuki- Trotter 
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formula we can express the fermionic partition function as 

Zf = Tr[e-^/^" 

= ('^il [1 -e^^i] \n'2){n2\--\n2d){n2d\ [1 - e^^2d] \n2d+i) 

ni,n2... 

('^2rf+i|[l - e^i] \n2d+2){n2d+2\---\n2Md){n2Md\ [1 - eifad] ìniXG) 

where the imaginary time extent 1/T has been divided into M equal steps 
of size e = 1/TM and each of these steps have been further divided into 2d 
time shces in each of which one of the ifj's act individually. We have also 
used the complete set of fermion occupation number states to evaluate the 
trace. 

The confìguration of fermion occupation numbers [n] = {nx,t},t = 
1,2, ...,2dM yields a fermion world line confìguration. The path integrai 
of the model is given by 

^/ = E SignM (7) 

[n] 



where the magnitude of the Boltzmann weight W[n] is the product of 
the magnitude of transfer matrix elements and the sign of the Boltzmann 
weight Sign [ri] is the product of their signs. Practically, if we ignore the 
anti-commutation relations between fermionic operators on different sites, 
the magnitude of the transfer matrix element and hence W[n\ does not 
change. On the other hand Sign[n] turns into a product of signs of the new 
transfer matrix elements times a global sign factor that takes into account 
the signs arising due to anti-commutation relations that were ignored in 
while calculating the transfer matrix elements. This global sign factor is 
topological in origin and can be found by tracking the permutation of the 



fermion world lines in time when the fermions are conserved |12|. The sign 



is positive for an even permutation and negative for an odd permutation. 



3 Meron Cluster Approach 

The world-line approach to fermionic path-integral cannot be used to de- 
sign algorithms because the Boltzmann weight is not positive definite and 
the correct probability density to generate world-line configurations is not 
known. One typically needs an exponentially large amount of statistics 
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to evaluate an expectation value using Monte-Carlo techniques. This is 
referred to as the sign problem. For example it is easy to check that the 
naive approach to evaluate 

{0) = ^^0[n]Sign[n]W[n] (8) 

^/ [n] 

where the partition function Zf is given by 

Zf ^Y^Sign[n]W[n] (9) 

[n] 

such that W[n] > and Sign[n] = ±1 suffcrs from a sign problem. 
Fortunately, in most physically intcrcsting problcms there is a lot of frcc- 
dom to choose the variables in which to express the partition function and 
evaluate the observables. Using this freedom some times one can be clever 
and find variables in which the Boltzmann weights turn out to be positive 
so that the sign problem is solved. 

Rccently, non-local cluster variables havc bccn succcssful is solving 
fermionic sign problems. The world-line partition function is first rewrit- 
ten by introducing new "bond" variables that connect lattice sites in ad- 
dition to the fermionic occupation variables that live on sites. Mathemat- 
ically this means 

Zf = J2 Sign[n] W[n] = ^ Sign[n, b] W[n, b] (10) 

[n] [n,b] 

where [n, b] refers to the new configuration of fermions and bonds. Bonds 
connect lattice sites into clustcrs and each bond configuration naturally 
divides ali sites of the lattice into a collection of clusters. New configu- 
rations can be obtained by reversing the fermion occupation on the sites 
associated with a single cluster. This is referred to as a cluster flip. 

Clearly, there is a lot of freedom in choosing Sign[n, b] and W[n, b] such 
that the partition function remains unchanged. However, if we restrict the 
choices such that 

(a) Cluster flips do not change W[n, b], i.e., W[n, b] = W[b], 

(b) Cluster fiips effect Sign[n, b] independently, 

(c) Cluster fiips can always produce a positive configuration, 
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then we can ensure a complete solution to the sign problem sirice it is 
possible to perforili an average over cluster flips which leads to 

Zf = T.^[b]mbl (11) 

[b] 

where 

J2 Sign[n, 6] = Siiìl[6] = 0, (12) 

clusterflips 

and Ne is the number of clusters. 

A cluster whose flip changes the sign of a configuration is called a 
meron. Properties (a) and (b) imply that such clusters identify two confìg- 
urations of equal weight and opposite signs and hence do not contribute to 
the partition function. On the other hand, meron clusters can contribute 
to observables. For example typically condensates get contribution from 
one meron sector, two point functions get contribution from zero, one and 
two meron sectors etc. An algorithm for such a problem must generate 
bond confìgurations with weight W[b] but suppress meron clusters. 



4 Spin-less Fermions 

Let US illustrate the ideas of the previous section using a simple example. 
Consider the nearest neighbor Hamilton operator introduced in section ^ 
with 

K;i = [-t(4c^+- + c+^fx) + + n^+i) +U{n^- 1)K.+ - - l; 

(13) 

with the constraint t = U/2 + fi. Here = clcx, is the fermion number 
operator. Figure |l| shows the non-zero transfer matrix elements [1 — eh^^y] 
in the occupation number configuration [n] as well as in the extended 
configuration of occupation numbers and bonds [n, b], such that partition 
function does not change. A typical configuration of fermion occupation 
numbers and bonds in one spatial and one temporal direction with M = 4 
is shown in Fig. ^ The shaded regions represent the interaction plaquettes 
each of which represent s a transfer matrix element. The global fermion 
permutation sign factor for this configuration is — 1. 

The signifìcance of the breakup chosen in Fig. |I| is that it yields a 
model that satisfies properties (a) and (b) of the meron cluster approach. 
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The first property is easily checked by noticing that the magnitude of the 
weights of the transfer matrix elements depend only on the bonds and not 
on the fermion occupation numbers. In order to check that the property 
(b) is also satisfied it is important to find the effect of a cluster fiip on 
Sign[n, 6]. It can be shown that if 

^side hops ^ , 

2 ~'~ cross hop over fiUed sites ~'~ temporal winding ^ ^ 

for a cluster is even only then the fermion permutation sign changes when 
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Figure 1: Interaction plaquettes and their weights that contribute to W[n] and 
W[n, b] so that Zf is preserved. It is assumed that t = U/2 + fi. 
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Figure 2: A typical configuration of fermions and bonds. 

that cluster is flipped. Details of the derivation of this formula can be 
found in []13|. The formula shows that the cross bonds can in principle 
induce dependence between clusters. However, the extra locai negative 
sign associated with the fully occupied cross bond in Fig. |I|, cancels such 
dependences. Finally for property (c) to be satisfied one has to either set 
/i = or [/ = 0. When /i = there are no cross bonds and one can 
flip ali configurations [n, b] to a configuration of staggered spatial fermion 
occupation which is static in time. This configuration is again guaranteed 
to be positive 0. On the other hand when U = there are no horizontal 
bonds and one can flip ali the configurations [n, b] to a configuration with 
ali sites are empty. This configuration is guaranteed to be positive H. 



5 Fermions with Spin 

Fermionic Hamiltonians that are of interest in both nuclear and condensed 
matter systems, involve internai degrees of freedom like spin and isospin. 
In order to illustrate how the meron cluster approach can be extended to 
such systems let us construct a model of fermions with spin which can 
be viewed as two layers of spin-less fermions. Instead of starting from 
the Hamilton operator and then constructing the cluster model as we did 
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spin "down" layer 





£a 

spin "up" layer 

Figure 3: Bond configurations and their weights of an interaction cube for a 
model with spin. 



in the case of spin-less fermions, let us construct the cluster model first. 
For simplicity, consider models where the bond configuration of Fig. ^ 
represents a typical configuration on each of the layers. Let the interac- 
tions between layers arise due to a constraints on how the bonds in the 
two layers are related. Fig. ^ shows the allowed bond configurations on 
each "interaction cube" of a simple interacting model and the magnitude 
of their weights. This means that every configuration is forced to contain 
identical clusters in the two layers. In addition, if the allowed fermion 
occupation numbers on the sites for a given bond configuration satisfy the 
rules of the previous section (i.e., the occupation numbers of sites con- 
nected by vertical and cross bonds must be the same and those connected 
by horizontal bonds must be opposite) and if ali these configurations have 
the same weight W[b] ( which is the product of the magnitudes of the 
interaction cube weights specified in Fig. |^) then the model automatically 
satisfies the property (a) of section ||. Since the clusters live on a single 
layer we can again use eq. (jl^) to determine the sign change due to a clus- 
ter fiip. In order for the cross bonds not to violate property (b), we again 
associate an extra locai negative sign with fuUy occupied cross bonds on 
each spin layer. It is easy to check that these restrictions automatically 
also satisfy property (c) and the cluster model thus becomes solvable with 
the meron cluster approach. 



As before the partition function of the model is given by eq. (11) with 
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Sign[6] given by eq. (|T2|). In the present case, the number of clusters Ne is 
always even with Nc/2 clusters in the spin up layer and an equal number 
of identical clusters in the spin down layer. When the quantity of eq. (|T^) 
is even for any of the clusters then Sign[6] = 0. Such clusters are the 
meron clusters. When there are no meron clusters in the confìguration 
then Sign[ò] = 2^'^ since each cluster can have two fermion occupation 
number configurations associated to it. Retracing the steps of the previous 
section in the reverse order, we can also determine the associated Hamilton 
operator. Remembering that 2dM is the number of time slices of the 
lattice and eM = 1/T, it can be shown that the model obtained in the 
limit e —>■ can be described by the Hamilton operator of eq. (|^) with 

[cz,T^S,T + Cj/,T^Ca;,t - n^, t - ny^^ + 1]} 
[Cx,l^Cy,l + Cy,i^c^,[ - 2K,i - 1/2) - 1/2) + 1/2] X 
[c.,T^c^,T + Cj;,T^c,,T - 2K,^ - l/2)(n^,T - 1/2) + 1/2]} . (15) 

Here the spin up and down fermions are created and annihilated by 
^x,]i^t,i Cx,],Cx,[ respectively and n^^] and n^^i refer to the corre- 
sponding number operators. 

The above Hubbard type model has a t/(l) fermion number symme- 
try. In three or more dimensions the breaking of this symmetry leads 
to superfluidity(or superconductivity when the symmetry is gauged). In 
two dimensions due to the Mermin Wagner theorem the symmetry can- 
not break spontaneously and the superfluid transition is driven by the 
Kosterlitz-Thouless phenomena |T^. It has been difìicult to study the 
predictions of universality in fermionic systems. This is the reason no 
known calculation exists that confìrms the predictions of the Kosterlitz- 
Thouless phenomena starting from a microscopie fermionic Hamiltonian. 
The above model for the first time provides an opportunity for studying 
such a superfluid transition using cluster algorithms. The simplest ob- 
servable relevant to this transition is the pair susceptibility which can be 
defined as 

^ = ±L r dt Tri e-(i/^-*)^ p+ e-*^ (16) 



ZV Jo 

with p+ = J2x '^I,t'^I,ì P^^^ creation and p" = (p^)^ the pair annihila- 
tion operators. In terms of cluster variables the susceptibility is propor- 
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Figure 4: Fair susceptibility ctS ci function of L. 



tional to the sum over the square of the size of certain clusters depending 
on the number of meron clusters in the configuration. The Kosterhtz- 
Thouless prediction says that 

^ °^ I Const. T > (1^) 

where rj{T) changes continuously from 0.25 at to at T = 0. Figure ^ 
shows the resuhs for the pair susceptibihty as a function of spatial size L 
in the above model with 6 = l,a = l,d = 2 and M = 20. As can be seen, 
T = 1.538, 1.429 are both above T^, although in the latter case lattices 
of size L = 128 are necessary to see the saturation. At T = 1.295 and 
T = 0.5 the power law fìts are extremely good over the entire range of L 
with powers 1.767(2) and 1.918(3). More details can be found in WU\. 



6 Model Extensions 

Although the meron cluster approach helps to find models without sign 
problems, it is also possible to find cluster models which seem not to 
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satisfy ali the properties of the meron cluster approach but stili do not 
suffer from sign problems. For example the model with the Hamilton 
operator H = Y.x,i hx,x+i + Ex hx, where 

hx = U (cx,^ - (^Cx,i - ^) + t^i^^,1 + "^^,1)- (18) 



will have the same partition function as eq. ( |Tl| ) except that Sign[ò] is no 
longer 2^^ lQ-^J^l jg replaced by 



(e/2d) U Se/ 4 



■Q Li^/2d) {-U 5c/4+Mfìc) _^ g(e/2d) {-U Sc/i-l^^c)'] ± 2e 

c e T layer 

(19) 

where Sq is the size of the cluster C and VLc is 2dM times the cluster's 
temporal winding number. The negative sign should be taken for a meron 
cluster. This factor is obviously positive for any if f/ < 0. Interestingly, 
it is also positive for /i < f//2 if 5 = and f/ > since then meron 
clusters always come in pairs. The proof also uses the fact that for ali 
clusters Sq > ^c- Thus we have found a repulsive Hubbard type model 
which does not suffer from a sign problem when formulated in the cluster 
approach for at least a limited range of chemical potentials. 



7 Conclusions 

We have sketched how a variety of fermionic partition functions can be 
written as a sum over cluster confìgurations with positive definite weights. 
The three properties of the meron cluster approach discussed in section ^ 
help find such models. However, it is possible to relax these properties and 
fìnd extensions in certain cases. One drawback of this approach is that the 
new models possess a complicated Hamilton operator. On the other hand, 
given the numerical efficiency of the algorithms that can be constructed 
for them it may stili be useful to study these models. For example this 
approach has led to the first precise confìrmation of universality arguments 
in fermionic systems [7-10]. Other applications in condensed matter and 
nuclear physics are being explored. 
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